1/f noise and avalanche scaling in plastic deformation 
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We study the intermittency and noise of dislocation systems undergoing shear deformation. Simula- 
tions of a simple two-dimensional discrete dislocation dynamics model indicate that the deformation 
rate exhibits a power spectrum scaling of the type l/f a . The noise exponent is far away from a 
Lorentzian, with a ~ 1.5. This result is directly related to the way the durations of avalanches of 
plastic deformation activity scale with their size. 
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I. INTRODUCTION 

Traditionally, the plastic deformation of solids has 
been thought of as a smooth process, so that the fluc- 
tuations can be neglected on large enough scales as the 
behavior becomes homogeneous. Contrary to this old 
paradigm, very recent evidence from experiments points 
out clearly that plastic deformation proceeds via bursts 
of dislocation activity [3, 0, 0, 0, S, B, % This activity, 
due to the long-range interactions, should have extended 
spatial and temporal correlations. Supporting this new 
picture, a number of simulation studies of simplified mod- 
els 0, H, 0, [l(| have demonstrated that such bursts or 
avalanches can be characterized by, apparently, scale-free 
size distributions. 

These advances bring plasticity close to a multitude 
of systems exhibiting "crackling noise" This is a 

generic property of many non-equilibrium systems, such 
as earthquakes [1 2 L Barkhausen noise in ferromagnets 
[IH, fracture [l4 Il5j. vortex avalanches in supercon- 
ductors [l6| and martensitic shape memory alloys [ItJ 
to name a few. The common feature among these are 
metastability and avalanches of activity that intercept 
quiescent periods. The realization that plastic deforma- 
tion has rich internal dynamics is of interest both for 
physics and materials science. Dislocations provide an 
interesting example of the physics of systems with long- 
range interactions competing with pinning from obstacles 
such as impurities, grain boundaries, and other disloca- 
tions [l|, [l8| ( "forest hardening" being a well-known ex- 
ample of the ensuing phenomena). 

Such avalanching systems are most simply studied by 
the one-dimensional activity time series V(t). In plas- 
ticity, V could be the global strain rate, or the acoustic 
emission activity during deformation, or in a constant 
strain rate experiment the shear stress. In this work we 
consider the temporal characterization of dislocation ac- 
tivity and show that its power spectrum (PS) exhibits 
what is known as 1//- noise, P(f) ~ 1// Q - This "flicker" 
noise is a well-known but relatively little understood phe- 
nomenon in various fields of science [191 ] . Then, we relate 
the scaling exponent a to the characteristic scaling of the 
avalanches. This has to be done with care due to the low- 



level background activity included in V. At the end, we 
are able to demonstrate that the noise of irreversible de- 
formation may follow generic, experimentally verifiable 
scaling laws. We use a simple two-dimensional (2d) dis- 
crete dislocation dynamics (DDD) model and study 
its behavior in the steady-state regime with a constant 
shear stress and a small (on the average constant) de- 
formation rate. While the statistics of avalanche sizes in 
this kind of models has been established to be of power 
law type @, Q, no temporal scaling analysis has been 
presented so far. 

The PS P(f) of a time series V(t) is the absolute square 
of the Fourier transform of V(t), 



P(f) 



J2nft 



V{t)dt 



(1) 



Since the time-time correlation function and the PS are 
related, the latter is a measure of temporal correlations 
in the system, a- values less than two indicate the pres- 
ence of complex time correlations. Recently it has been 
realized that under certain fairly general conditions the 
exponent a can be derived from the scaling exponents 
characterizing the avalanche distributions, the examples 
ranging from models of Barkhausen noise [2(J, to self- 
organized criticality [21] to fluid invasion into disordered 
media 22]. The intermittency is easier to understand in 
toy models as sandpile models of self-organized critical- 
ity [23[, where V(t) drops to zero between avalanches. 
In most experiments and more realistic model systems 
- such as here - there is a background (from noise, or 
processes that co-exist with the avalanches). Also, finite 
drive rates can lead to merging of avalanches, thus to 
problems in characterizing the underlying activity. 

The structure of the paper is as follows: In Sec. II 
we introduce the 2d DDD model, and in Sec. Ill the 
scaling relation relating the power spectrum exponent 
to the avalanche statistics is presented. Results of the 
numerical simulations are presented in Sec. IV. Section 
V finishes the paper with conclusions. 
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II. DDD MODEL 

Despite its simplicity, the DDD model we study here 
has been shown to capture many features of realistic plas- 
ticity, such as avalanches with scale-free size distribu- 
tions [2, .&l as well as the Andrade creep law |Ti . It is 
a two-dimensional model representing a cross section of 
a crystal with single-slip geometry. The dislocations are 
assumed to be straight parallel edge dislocations, with 
the dislocation lines oriented along the z axis. They 
glide along directions parallel to their Burgers vectors 
b = ±bu x , with u x the unit vector in the x-direction. 
It is thus sufficient to consider a two-dimensional system 
(i.e. the xy plane) with N pointlike edge dislocations 
gliding in the x-direction. Equal numbers of dislocations 
with positive and negative Burgers vectors are assumed. 
For simplicity other processes contributing to dislocation 
motion such as dislocation climb are not considered. 

An essential feature of the model is that the disloca- 
tions interact with each other through their anisotropic 
long-range stress fields 



a s (r) = Db 



x(x 2 — y 2 ) 



(x 2 + y 2 ) 



2\2 ' 



(2) 



where D = n/2w(l — v), with /i the shear modulus and v 
the Poisson ratio of the material. We assume overdamped 
dynamics with a linear force- velocity relationship, giving 
rise to the equations of motion of the form 



Xd Vr, 



s„b 



^ s m a s (r nm ) + a 



(3) 



where v n is the velocity of the n'th dislocation, Xd is the 
dislocation mobility, s n refers to the sign of the Burgers 
vector of the n'th dislocation and a is the external shear 
stress experienced by the dislocation. These equations 
are put into a dimensionless form by measuring lengths in 
units of b, times in units of l/(x<iDb) and stresses in units 
of D. Due to the long range nature of the dislocation- 
dislocation interaction, periodic boundary conditions are 
implemented with an infinite sum of images [25) . The 
positions r„, n = 1,...,N, of the N dislocations as a 
function of time are computed by integrating the equa- 
tions of motion (|3|) numerically by using an adaptive step 
size fifth-order Runge-Kutta algorithm. 

For small distances, Eq. @ ceases to be valid. Thus 
for distances smaller than 2b from the dislocation, we 
set the dislocation stress field to zero. This procedure 
removes the unphysical singularity in Eq. j2|). Fur- 
thermore, when the distance r nm between two disloca- 
tions with Burgers vectors of opposite sign gets small, 
i.e. r nm < 2b, we employ a phenomenological annihila- 
tion reaction by simply removing them from the system. 
To compensate, and to include dislocation multiplication 
(as in real plasticity through Frank-Read sources), we 



introduce a mechanism to create new dislocations. The 
system is split into smaller cells and then one monitors 
both the local stress and the number of pinned disloca- 
tions (i.e. those moving slower than a threshold velocity) 
in each cell. Given pinned dislocations and that the local 
stress exceeds a threshold value, a dislocation pair with 
opposite Burgers vectors is generated with a probabil- 
ity proportional to the absolute local stress value. This 
means that on the average it takes a finite time for the 
source to create a new dislocation pair, as is the case 
with real dislocation sources. The two new dislocations 
of opposite sign are inserted into random locations inside 
the neighboring cells of the one containing the activated 
source. This is done with the constraint that their com- 
bined stress field must decrease the magnitude of the lo- 
cal stress at the source location, while elsewhere in the 
system further dislocation activity may be triggered, re- 
sulting in a correlated sequence of dislocation activity. 



III. SCALING OF THE POWER SPECTRUM 

Consider a bursty time series V(t) consisting of tem- 
porally separated avalanches. The usual definition of an 
avalanche is a connected sequence of values of V(t) ex- 
ceeding some threshold value V t h, to e.g. subtract uncor- 
rected background noise. If an avalanche starts at t = 
and ends at t = T, the size s of an avalanche of duration 
T is defined as s(T) = ^[V{t) - V th \dt. Assume that 
the average avalanche size (s(T)) of avalanches of a given 
duration T scales as a power law of the duration, 



(s(T)} - T 7s 



(4) 



and that the avalanches are self-similar so that the aver- 
aged avalanche shapes V(T,t) of avalanches of different 
durations T can be collapsed onto a single curve by using 
the ansatz 



V(T,t) = T^- 1 f shape {t/T) 



(5) 



Here, f s hape{x) is a scaling function. The total energy is 
obtained as the 9 — component of the stationary time- 
time correlation function, defined by C{9) — J V(t)V(t + 
0)dt. By cosine transforming the time-time correlation 
function C(9\s) of avalanches of a given size s, one obtains 
the scaling form 



(6) 



for the corresponding energy spectrum. The scaling 
of the total power spectrum then follows by averaging 
E(f\s) over the avalanche size probability distribution 
D s (s), e.g. a power law D s (s) ~ s~ Ts cut off at s = s*, 
so that 



p(f) = r lBt{ *- TB) 



dxx 



2-t s 



9e{x). (7) 
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The result will depend on the value of r s as well as on the 
form of the scaling function Qe(x) [IS HH- In the case 
g E {x) ~ 1/sg @, El and T s < 2 (for the DDD models 
at hand, the latter condition seems to be fulfilled, with 
t s 1.6 8]) the power spectrum scales as 

P(/) ~ (8) 

This indicates a scaling relation a = 7„t, which possibil- 
ity we next check with simulations. 

IV. RESULTS 

A natural choice for the time series is the "collective 
velocity" of the dislocations, V(t) — J2i \ v i\- It is propor- 
tional to the energy dissipated per unit time by the dis- 
location system. Thus the statistical properties of V(t) 
could be related to the acoustic emission (AE) statistics, 
as suggested by a number of authors, see Ref. An- 
other possibility is V s (t) = J2i ^i v ii i- e - the global strain 
rate. We mainly present results here for the first choice. 
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FIG. 1: An example of the time series V(t) = £V \vi\ from 
a system of linear size L = 3006. a = 0.03. Time is re- 
ported in units of l/(xdDb) and V in units of XdDb 2 . In the 
steady state, the system contains on the average about 600 
dislocations. 

The numerical simulations are started from a ran- 
dom initial configuration of No dislocations with equal 
numbers of positive and negative Burgers vectors (here 
L = 2006 and L = 3006, with A = 600 and A = 1000, 
respectively). The system is first let to relax in the ab- 
sence of external stress, until it reaches a metastable ar- 
rangement with N < Nq. Then a small constant ex- 
ternal stress is applied and the evolution of the system 
is monitored. In the absence of dislocation multiplica- 
tion, the system undergoes a "jammin g tr ansition" at a 
critical value a c of the external stress [24|. For a = a Cl 
the strain rate decays as a power law in time, with the 
Andrade's power law creep exponent close to —2/3. For 
a > <j c w 0.01, the system crosses over to a constant 
strain rate regime (linear creep) which we study next. 

With the dislocation multiplication turned on, the sys- 
tem displays intermittent avalanche-like bursts of dis- 
location activity. The number of dislocations fluctu- 
ates around N « 300 for L = 2006 and N w 600 for 



L = 3006. Fig. [T] shows that V(t) displays intermittent 
behavior, consisting of avalanches. These are defined as 
a connected sequence of V(t)- values exceeding a thresh- 
old value Vth and have a wide range of sizes. Here, we 
monitor V(t) for external stress values close to but above 
the critical value o~ c . 




FIG. 2: (Color online) The average scaled avalanche shape for 
different threshold values Vth f° r a = 0.03 and L = 3006. For 
small threshold values the avalanches are asymmetric with the 
asymmetry decreasing with Vth- The thick lines with circu- 
lar symbols correspond to an average over of the order of 10 3 
avalanches, which have been scaled according to Eq. ((5} be- 
fore averaging. The thin solid lines correspond to avalanches 
of different limited duration ranges. The inset displays the av- 
erage creation and annihilation rates (r c and r a , respectively) 
of dislocations during an avalanche, with Vth = 5. 



Fig. [2] exhibits the average avalanche shape (one needs 
a minimum duration, of the order of 10, for scaling), 
for varying Vth- They appear to agree with Eq. ([5]). 
The data has been obtained by scaling all individual 
avalanches in a given duration range according to Eq. 
([5]). The threshold value Vth does not have any effect 
on the s(T) statistics (given that Vth is high enough such 
that a "noise level" corresponding to incoherent motion of 
dislocations is not taken into account), but the avalanche 
shapes appear to depend on Vth- For small Vth -values, 
avalanches are clearly asymmetric to the left (in agree- 
ment with experiments 26]), but become more symmet- 
ric upon increasing Vth ■ This asymmetry is also manifest 
in the creation rate of dislocations during an avalanche 
(inset of Fig. [5]). In Fig. [3] we show the energy spectra 
computed using the Lomb periodogram, scaled according 
to Eq. ©. The behavior of the scaling function gs(x) 
indicates that Qb{x) ~ X/x. 

The main result is shown in Fig. [U where we compare 
the scaling of the total power spectrum, from V(t), with 
that of (s(T)). We observe a PS of the form P(f) ~ f~^ Bt 
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FIG. 3: (Color online) The scaled energy spectrum of dislo- 
cation avalanches, showing the form of the scaling function 
gE(x), which is consistent with the <?b(x) ~ 1/x behavior, 
indicated by the thick solid line. The data has been obtained 
by computing the energy spectrum E(f\s) of each avalanche 
for a = 0.03 and L = 3006, normalizing by s 2 , and averaging 
over all avalanches longer than a threshold duration. 



spanning almost two decades, with "f s t ~ 1.5. Thus the 
PS of dislocation activity is related here to the intrinsic 
scaling of the avalanches. The extension of the scaling 
regime increases with L, the system size (Fig. [5]), with 
a cut-off frequency f C ut-of / roughly equal to the system 
size dependent inverse duration of the longest avalanche, 
f cut-off ~ T^ X {L). The absence of scaling for the very 
highest frequencies (/ > 0.1) is due to a finite crossover 
time, before the avalanches have a self-similar structure. 
Similarly to the exponent of the avalanche size distribu- 
tion 0, the exponent of the power spectrum is interest- 
ingly insensitive to the value of the external stress. In the 
inset of Fig. 01 we consider the PS of the total strain rate. 
While the scaling region appears somewhat narrower in 
this case, the exponent a is observed to be unchanged 
from 1.5. 

These results are largely independent of the details of 
the model, such as the threshold value for the local stress 
to create new dislocations. We have checked that the 
same results can be recovered even in the case with no 
dislocation multiplication, where the avalanches are due 
to different threshold mechanisms, such as unpinning of 
dislocation dipoles. As this approach suffers from the fact 
that simulations are computationally more demanding 
due to the longer time scales one must reach, we have here 
restricted ourselves to considering only the case where 
dislocation multiplication does occur. 



Knowing that the the avalanches are self-affine one can 
also compute "f s t from the other known exponents. Con- 
sider the probability Dy{V\s) that a value V will occur at 




f, 1/T 

FIG. 4: (Color online) A comparison of the scaling of the 
power spectra (solid lines) with that of (s(l/T)} (dashed 
lines), for a system of linear size L = 3006 and various a- 
values. The thick solid line corresponds to a — 7 3t = 1.5. 
The inset shows the power spectrum of the strain rate, for 
a = 0.04. The results are averaged over several random ini- 
tial configurations. Both / and 1/T are reported in units of 
XdDb. 

some point during an avalanche of size s. Assuming scal- 
ing and requiring normalization, Dy(V\s)dV = 1, as 
well as (V) = s/T ~ s 1 " 1 ^.., 

one obtains (see also I20H ) 

D v (V\s) = V-'MVs 1 ^- 1 ), (9) 

where fv(%) is a scaling function. The total probability 
Dv(V) is then obtained by integrating Dv(V\s) over the 
avalanche size distribution D s (s) ~ s ~ Ts , giving rise to 

D V (V) ~ = V~ TV . Miguel et al. @ studied 

the statistics of the quantity E — (J^ |wi|) 2 , and found 
that its distribution behaves like De(E) ~ E~ TE , with 
te ~ 1.8. From te = \(jv + 1) one ma y solve for j st , 
resulting in (with t s « 1.6 Q) 



in reasonable agreement with our results above. 

V. CONCLUSIONS 

In this work, we have demonstrated that a simple two- 
dimensional dislocation system exhibits 1 / /"-noise, so 
that the relation a = j s t can be explained by the scal- 
ing properties of the avalanche dynamics. We have also 
considered the PS of the global strain rate, with similar 
conclusions. These findings lend themselves to experi- 
mental tests, e.g. as in the creep experiments of ice single 
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FIG. 5: (Color online) Power spectra for two different sys- 
tem sizes, L = 200b and L — 300b, corresponding to on the 
average 300 and 600 dislocations in the steady state. The 
extension of the scaling region of the PS is seen to increase 
with L. The results are averaged over several random initial 
configurations. The frequency / is reported in units of XdDb. 

crystals in experiments on deforming metallic single 
(micro)crystals and could also be considered in ex- 

periments on colloidal crystals [28|. In AE experiments 
performed with three-dimensional systems, values consis- 
tent with mean field exponents have been reported 0,[23] , 
suggesting a — "f st — 2 (see Eq. (fTTJ)) ). It is thus possi- 
ble that to get quantitative agreement with experiments, 



the full three-dimensional problem should be studied, e.g. 
by means of simulations of a three-dimensional disloca- 
tion dynamics model. Other issues not considered in the 
present study include the possible relevance of quenched 
disorder e.g. in the form of forest dislocations, which 
would provide strong pinning centers to resist dislocation 
motion. One should note, however, that also the present 
model contains a pinning force landscape generated by 
the dislocations themselves 0). 

In polycrystals, avalanches interact with the grain 
boundaries, which will probably lead to a size-dependent 
avalanche shape (29[. This means that Eq. is no 
longer directly applicable, presenting an interesting the- 
oretical question. Our theory assumes negligible corre- 
lations between avalanches. In experiments, dislocation 
avalanches have been found to exhibit a tendency to clus- 
ter in time, such that a "mainshock" is typically followed 
by a sequence of few "aftershocks" [3, While this 
may modify the low frequency part of the PS, the high 
frequency part, corresponding to correlations within in- 
dividual avalanches, is still expected to scale according 
toEq. ®. 

Finally we note that in addition to materials with crys- 
talline structure, one field where the use of spectral tools 
and the study of avalanches might be used to charac- 
terize the spatio-temporal behavior is the plasticity of 
non-crystalline media, e.g. amorphous glasses, where the 
localization and intcrmittency of plastic events has been 
recently demonstrated in simulations [3(| HH . 
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